arXiv:1506.06758v2 [cond-mat.supr-con] 24Jun2015 


Terahertz time-domain spectroscopy of transient metallic and superconducting states 

J. Orenstein^’^ ’0 and J. S. Dodge^ '^’S 

^Department of Physics, University of California, Berkeley CA 94720, USA 
^Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley CA 94720, USA 
^Department of Physics, Simon Fraser University, Burnaby, British Columbia, V5A 1S6, Canada 
^Canadian Institute for Advanced Research, Toronto, Ontario, M5G 1Z8, Canada 

Time-resolved terahertz time-domain spectroscopy (THz-TDS) is an ideal tool for probing photoinduced 
nonequilibrium metallic and superconducting states. Here, we focus on the interpretation of the two-dimensional 
response function that it measures, examining whether it provides an accurate snapshot of the instan¬ 

taneous optical conductivity, o{(0\t). For the Drude model with a time-dependent carrier density, we show 
that L{(0\t) is not simply related to G{(0\t). The difference in the two response functions is most pronounced 
when the momentum relaxation rate of photocarriers is long, as would be the case in a system that becomes 
superconducting following pulsed photoexcitation. From the analysis of our model, we identify signatures of 
photoinduced superconductivity that could be seen by time-resolved THz-TDS. 


PACS numbers: 78.47.-p, 78.56.-a, 42.65.-k, 74.25.N- 


I. INTRODUCTION 

Time-domain terahertz spectroscopy (THz-TDS) probes 
the optical conductivity of metals and superconductors by 
measuring the current transient induced by a subpicosecond 
electric field pulse. Time-resolved THz-TDS exploits the 
short duration of the THz probe to detect how the conductiv¬ 
ity changes in response to pulsed photoexcitation after a con¬ 
trolled delay—providing a way to take snapshots of the optical 
conductivity with picosecond time resolution. This scheme 
has been used effectively to measure transient photoconduc¬ 
tivity in a wide variety of bulk and nanostructured semicon¬ 
ductor systems^i and its range of application is growing^^Tii 
Recently, THz-TDS spectra of the photoexcited normal state 
of high-Tc cuprates and K 3 C 60 have been presented as evi¬ 
dence of transient photoinduced superconductivity, because 
they resemble the equilibrium conductivity spectra obtained 
well below the superconducting transition temperature 
However, the straightforward interpretation of THz-TDS as 
a snapshot of the optical conductivity spectrum breaks down 
when the characteristic relaxation times are comparable to the 
photoexcitation delay, because the Fourier transform involves 
times that precede photoexcitation 

In order to assess the importance of this issue for in¬ 
terpreting measurements on transient metallic and super¬ 
conducting states, we analyze a simple model of tran¬ 
sient photoconductivity^i^-— For photoexcitation at r = 0, 
we compare the time-dependent instantaneous conductivity 
change, 5a{(0\t), with the response function Tl{(0\t) mea¬ 
sured in time-resolved THz-TDS. We find that Z shows large- 
amplitude deviations from 5(7 when the photocarrier Drude 
scattering time is longer than t^di regime that may be acces¬ 
sible in a transient photoinduced superconductor. Our analy¬ 
sis offers improved guidance on how and when time-resolved 
THz-TDS spectra can be interpreted as a conductivity snap¬ 
shot of a transient state. 

We consider a pump-probe experiment in which a material 
is photoexcited by a strong pump pulse with intensity profile 
I{t) « /o5(r), and the current induced by a THz-frequency 
probe field is measured at time t. The photoinduced change 


in the current is then 

5J{t) = Jj E{t-x)I{t-Xe)a^^\x,Xe)dTdXe, ( 1 ) 

where E{t — t) is the electric field of the THz probe at the 
moment T before the observation time, and is a 

third-order susceptibility, in sense that its associated current 
is proportional to both the THz probe field and the pump in¬ 
tensity (two powers of field). By defining 

/ CX3 

l{t-Xe)(y^^\x,Xe)dXe, ( 2 ) 

-oo 

we can rewrite Eq. (ID as 

5Jit)= [ E{t-x)L{x,t)dx, (3) 

which has the usual linear response form—except, crucially, 
that the two-dimensional response function Z(T,r) has an im¬ 
plicit dependence on I{t) and lacks time invarance.-“— Kindt 
and Schmuttenmaer (KS) pointed out that could be 

readily obtained with time-resolved THz-TDS by indepen¬ 
dently controlling the delay of the THz probe field with re¬ 
spect to photoexcitation, as illustrated in Fig. [T] and described 
more thoroughly in Sec.HIIU^ 

With a measurement of Z(T,r), it is straighforward to com¬ 
pute the Fourier transform Z((0;t). The question that arises 
then is the following: whether, or in what limit, can Z(ci>;t) 
be considered equivalent to the pump-induced change in the 
instantaneous optical conductivity, Sa((0;t)?^~— Indeed, 
whether a response function Sa(co;t) exists at all for an ar¬ 
bitrary nonequilibrium system is itself a problem, one we do 
not consider here. Instead, we consider a specific example of 
a class of optically pumped systems for which the concept of 
an instantaneous linear response function should be applica¬ 
ble. Following photoexcitation, this class of systems passes 
through a continuous sequence of quasi-equilibrium states, 
in which the distribution functions of electrons, phonons, 
magnons, etc. can be described by quasi temperatures and 
chemical potentials. For such systems it should be possi¬ 
ble to define the response function Sa(co;t) as the change in 
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FIG. 1. (Color online) (a) Schematic illustration of the pulse se¬ 
quence and time delays in a time-resolved THz-TDS measurement 
performed on the model system described in the text. The red lines 
represent the photoexcited carrier density that is generated by a pump 
pulse at t = 0, then decays on a timescale l/F = 20 ps. The black 
lines represent the electric field of the THz probe pulse, shown for 
increasing time delay, tp (marked by x), in the three panels. Blue 
lines depict the corresponding nonequilibrium current, and blue dots 
mark its value at = 2 ps. The dotted line shows the impulse re¬ 
sponse function a{ts —t) for a Drude metal. The truncation of the 
impulse response at t = 0 differentiates from 5a{T,ts). (b) 

Nonequilibrum current as a function of tg — tp for fixed tg =2 ps and 
1/r = 20 ps. With increasing momentum relaxation time, Tg, a pulse 
of nonequilibrium current centered on ^^ = 0 grows in amplitude. 


the equilibrium a{co) that would be measured if the quasi¬ 
equilibrium state at time t were entirely metastable. Below 
we show by example that although 5a{cou) can be well de¬ 
fined, it is not, in general, equivalent to the function ^{cou) 
measured by time-resolved THz-TDS. 


II. EXAMPLE OF PHOTOCARRIERS WITH DRUDE 
RESPONSE 


We treat the Drude model for nonequilibrium photocarriers, 
but it is helpful to start with elementary equilibrium relation¬ 
ships. For a system of n carriers with Drude scattering rate 7 , 
the current is given by 


/ oo 

E{t -x)a{x)dT, 

-00 

with impulse response 


(4) 


<y{x) = —©(T)exp(-7T). (5) 

m 

The Drude conductivity spectrum is then just the Fourier 
transform of Eq. ([5]): 


g{(o) 


ne^ 1 
m y—i(0 


( 6 ) 


We turn now to the nonequilibrium case. If photoexcitation 
creates Sn(0) carriers that recombine at a rate F, then the in¬ 
stantaneous conductivity is given by the Drude formula with 
n(t) = 5n(O)0(r)exp(—Fr), 


S(j(co;t) = exp(—rt) 
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y-ico 


(7) 


or, in the time domain, 


5a{T,t) = 0(r) 


m 


exp(—F r)0(T) exp(—7 t). 


( 8 ) 


We emphasize here that in both 5a{T^t) and 5cr(()L);r), the de¬ 
pendence on t is entirely through the state variable n{t), which 
is then assumed constant when we consider the dependence of 
5(7 on the dynamical variables T and (0. 

By contrast, E(T,r) has an explicit dynamical dependence 
on t. Because of the relative simplicity of our model, we can 
determine E(t, r) directly by integrating the classical equation 
of motion for the photocarriers. As shown in Appendix!^ we 
obtaini^ 


Z(T,r) = 0(^ — T)0(r) exp(—Fr)0(T)exp(—7 t). 

m 

(9) 

Eq. ® is identical to Eq. ([8]) except for the additional Heav¬ 
iside function, 0(r — t), which prevents carriers from con¬ 
tributing to the integrand in Eq. (O before they are created at 
r = 0 (or T = t). 

Substituting this form for E(T,r) into Eq. ^ yields. 


8J{t) = 0(f) exp(-rf) 

m 

X f E{t - x)&{t - x)&{x)txp{-Yx). (10) 

J -00 

In the following section we show how time-resolved THz- 
TDS is applied to measure E(T,r), and compare 'L{cou) with 
5G{co;t). 
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III. TIME-RESOLVED TERAHERTZ SPECTROSCOPY 

As a preliminary step we review the methodology of THz- 
TDS as applied to the equilibrium optical conductivityTHz- 
TDS effectively measures the current at a time ts that is in¬ 
duced by a THz electric field pulse, E{t) = Ep{t — tp), cen¬ 
tered on a time tp’. 

/ oo 

^pij's tp ( 11 ) 

The induced current J{tsdp) is inferred from a time-resolved 
measurement of the electric field reflected from, or transmit¬ 
ted through, a medium under test. It depends only on the dif¬ 
ference ts — tp between the sampling and probe arrival times, 
so it can be measured by scanning either 4 or tp. Using the 
convolution theorem, a{co) can be obtained by Fourier trans¬ 
forming J{tsdp) along the Ar = 4 — tp direction. 

Time-resolved THz-TDS focuses on the change in the re¬ 
sponse to the THz probe that is induced by a pump pulse. 
The sequencing of pump and THz probe pulses is illustrated 
in Fig. [Ua). Using Eq. (fTOl) . the nonequilibrium current for 
Drude photocarriers is 

8J{ts,tp) = exp{-rts) 

m 

X f Ep{ts-lp-t)&{ts-t)&{t)e\p{-Yt)dT. (12) 

J -00 

Note that 5f{tsdp) is now a function of both the sampling 
and probe arrival times, rather than just their difference, as 
in equilibrium THz-TDS; this reflects the breaking of time- 
invariance by the pump. An example of the transient nonequi¬ 
librium current calculated from Eq. (O is illustrated in 
Fig-Eb), which shows 5J{tsdp) as a function of ts — tp with 
ts = 2 ps, for several values of = I/ 7 . For = 0.3 ps, the 
current follows the electric field pulse, with a lag and slight 
distortion caused by convolution with the Drude impulse re¬ 
sponse. As Ts increases, a feature emerges near = 0 that has 
a completely different origin, illustrated in the tp = {) panel of 
Fig-Ea): here, only part of the THz field can induce current, 
creating an unbalanced current impulse that persists until the 
measurement time ts. When Ts as in a superconductor, 
these two features have equal and opposite magnitude; how¬ 
ever, their separation in time will vary with the observation 
time ts. 

The procedure introduced by KS is to scan tp at fixed ts, 
such that the integral expression for the current retains the 
form of Eq. (fTTI) . but with the two-dimensional nonequilib¬ 
rium response function in place of g{t), as in Eq. (E- 

If we now Fourier transform Z(t, 4) with respect to T, we get 

E(a);4) = exp(-r?^)—U 

m y-ico 

X [1 -exp(-74)exp(/a)4)], (13) 

which differs from the instantaneous conductivity 5(7(to; 4 ) in 
Eq. (E by the term in square brackets. This result is consistent 


with earlier results on the nonequilibrium Drude modelri“— 
expressed in a way that allows more immediate comparison 
with experiment. 

We have assumed impulsive excitation, I{t) ^ IoS{t), that 
causes the conductivity to change abruptly. When the excita¬ 
tion pulse width T^ cannot be neglected, we expect the oscilla¬ 
tions in E{cO’Js) to become damped as COT^ ^ 1. Another fac¬ 
tor that can lead to damping occurs if the photoinduced change 
in conductivity has a finite risetime. This would be reflected 
in the dependence of the third-order susceptibility Te) 

on Te in Eq. (E- In the presence of either form of broadening 
of the step-function change in conductivity, the deviation of Z 
from G will be most pronounced at low frequencies, which is 
indeed observed^2Jl 


IV. DISCUSSION 

Eq. ([T^ greatly clarifies the conditions under which the 
time-resolved THz-TDS spectrum Z((0;4) approximates the 
instantaneous linear response, SG(co;ts). First, it shows that 
photocarrier recombination simply rescales the overall spec¬ 
trum by exp(—rts), so the measurement fidelity is not fun¬ 
damentally limited by the recombination time. The critical 
parameter is yts, the product of the momentum relaxation rate 
and the sampling time. 

The crossover in the nature of the spectra at yts ^ 1 is il¬ 
lustrated Fig.Ea), in which we plot Z(co;ts) as a function of 
0) for several values of Ts, with 4 fixed at 2 ps. The instan¬ 
taneous Drude conductivity SG(co;ts) for the same values of 
Ts is shown as dotted lines for comparison, and spectra with 
different Ts are normalized to SG(co;ts =0) to illustrate the 
variation in frequency dependence. 

When ts/Ts is large, the component of the current sampled 
at ts from carriers accelerated at r = 0 is exponentially small. 
Consequently, Z(co;ts) asymptotically approaches the instan¬ 
taneous conductivity, which has the Drude form in our simple 
example. On the other hand, if Ts is comparable to or greater 
than ts, a component of the current that would be present if the 
state were metastable is cutoff at t = 0. The temporal cutoff 
generates oscillations with period 271/ts along the frequency 
axis of Z((0;r^), which is clearly no longer simply related to 
SG(co;ts). 

Fig. E shows this behavior in more detail as it would appear 
for photocarriers with = 1 ps and l/E = 0.5 ps. The imag¬ 
inary part of E(cD;ts) has oscillations that appear as ridges 
along constant values of cots, with an amplitude that decays 
exponentially with 74 . Similar oscillations were observed 
previously in both measurementsi^ and simulationsand 
Eq. [ 13 ] clarifies their origin. 

Finally, we consider what might be observed in a time- 
resolved THz-TDS measurement in which a transient super¬ 
conducting phase is generated at r = 0 by a laser pulse. De¬ 
veloping a phenomenological description of transient super¬ 
conductivity is clearly not as straightforward as modeling 
a transient photoconductor. One approach that is directly 
amenable to our analysis is based on the two-fiuid model, 
which describes the current response in terms of normal fluid 
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FIG. 2. (Color online) (a) Real and (b) imaginary part of E(ft);4) 
as a function of ft) for several values of with 4 fixed at 1 ps and 
r = 0. The Drude conductivity for same values of is shown as dot¬ 
ted lines for comparison. Spectra with different Zg are normalized to 
Oo = n{0)Tse^Im to illustrate the variation in frequency dependence. 
When tg/Zg is large Z{(o;tg) asymptotically approaches the instanta¬ 
neous (Drude) conductivity. For Zg > tg the temporal cutoff at ^ = 0 
generates oscillations, with period 2%/tg, along the frequency axis of 

^{(0\tg). 

and superfluid components, with spectral weights Une^ /m and 
UgC^ /m respectively. The normal fluid conductivity is de¬ 
scribed by the Drude spectrum, while the superfluid compo¬ 
nent is characterized by an inflnite momentum relaxation time. 
We can then formulate photoinduced superconductivity as the 
generation of superfluid spectral weight /m by transfer 
from the normal fluid. If we assume that the photoinduced 
superfluid has a lifetime l/F, then 

Ar? ^ 
m 

X |2[1 _exp(!(0/i)] - 2 [1 _exp(/wfi)]|, (14) 

with a = (0 + ij- Fig. m illustrates the spectra predicted by 
Eq.m for several values of tg, with =I ps and Zr =0.5 ps. 
The spectral shape is dominated by underdamped oscillations 
that originate from the sharp cutoff in the time-domain re¬ 
sponse at r = 0, as discussed above. We note that in this de¬ 
scription of transient superfluidity, Z(ft);^^) never approaches 
5(7 (ft); 4 ) because of the undamped contribution to f{tg^tp) 
from super carriers generated aX t = 0. With an overall sign 
change, Eq. ([TO should also describe Z(ft);4) when photoex¬ 
citation suppresses superconductivity, for example by trans¬ 
ferring spectral weight from the superconducting carriers to 
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FIG. 3. (Color online) Simulation of time-resolved THz-TDS re¬ 
sponse function for photocarriers whose dynamics are characterized 
by a time-independent scattering rate, (a) Imaginary part of Z{(o;tg), 
with = 1 ps and l/F = 0.5 ps, in both 3D and contour plot rep¬ 
resentations. Contour levels indicate 0.01 steps from 0.01 to 0.1, in¬ 
clusive. (b) Spectra of 'L2{<x>’Js) (solid lines) for several values of tg, 
with 5G2{co;tg) (dotted lines) shown for comparison. The spectrum 
for each tg is denoted by color in the legend; both Z 2 and <72 decrease 
with increasing tg, as the carriers decay. All spectra are normalized 
to (7o = dn{0)e^Zg/m. 


quasiparticles above the gap^^ 

As mentioned above, different phenomenological descrip¬ 
tions of a transient superconducting state are possible. For ex¬ 
ample, rather than generating a fully coherent superfluid com¬ 
ponent, one could imagine that the effect of the pump pulse 
is to suddenly increase the momentum relaxation time of the 
entire electron fluid to some large but still flnite value. The 
system would subsequently return to equilibrium through the 
decay of partial coherence and recovery of the normal state 
Zg. We describe the TD-THz response for the case of a time- 
varying momentum relaxation rate in the Appendix, where we 
obtain the response function E(T,r) in terms of Y{t). Unlike 
the models considered above, we believe that the instanta- 
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FIG. 4. (Color online) Simulation of time-resolved THz-TDS re¬ 
sponse function for a model of transient photoinduced superconduc¬ 
tivity. (a) Imaginary part of Z{co;ts), with = 1 ps and l/F = 0.5 ps, 
in both 3D and contour plot representations. Contour levels indicate 
0.005 steps from -0.01 to 0.04, inclusive, (b) Spectra of L2{co',ts) 
(solid lines) for several values of 4, with 5G2{co;ts) (dotted lines) 
shown for comparison. The spectrum for each tg is denoted by color 
in the legend; with increasing tg, G 2 decreases uniformly as the car¬ 
riers decay, while Z 2 shows oscillations about Z 2 ~ 0, with period 
In/tg and an amplitude that decays with both co and tg. All spectra 
are normalized to Aoo = l^ne^Xglm. 


neous conductivity, a{cou), is not well-defined in a system 
where / depends explicitly on the time. Still, ^^{cou) remains 
a valid response function, and can exhibit features that are 
similar to those predicted for the two-fiuid model. In this case, 
the deviations from a Drude spectrum will smaller than in the 
fully coherent superconductor, particularly if the maximum 
tg reached by the partially coherent state does not exceed the 
sampling time. 


V. SUMMARY AND CONCLUSIONS 


Time-domain terahertz spectroscopy provides an elegant 
method for doing time-resolved photoconductivity measure¬ 
ments. In this paper we focused on whether the response 
function ^^{cO’Jg) that is typically measured in time-resolved 
THz-TDS can be interpreted as the photoinduced change 
5a{co;tg) in the instantaneous optical conductivity. Within 
a simple model, we showed that £(( 0 ;^^), is never equivalent 
to 5 cr(m;r 5 ), although we also found that ^{cO'Js) approaches 
S(j(co;tg) asymptotically in the limit yfg 00 , where / is the 
Drude relaxation rate of the nonequilibrium carriers. In this 
limit, the current measured at fg has an exponentially small de¬ 
pendence on the field applied at t <0; that is, before the pump 
pulse arrives. In the opposite regime, ytg < 1, we presented 
an analytic expression that shows that Z(cD;tg) and S(j(cD;tg) 
are entirely distinct response functions for the nonequilibrium 
Drude model. Here, the absence of nonequilibrium carriers 
for r < 0 creates a current imbalance for fields applied near 
t = 0 that persists until the measurement time. Neverthe¬ 
less, we believe that even in this regime, information about 
the number density and mobility of the photocarriers can be 
obtained by comparing L(co;tg) with theoretical models^"— 

The effort to better understand the relationship of E(co;tg) 
to S(7(co;tg) was largely motivated by experiments reporting 
photoinduced transient superconductivity in cuprates and in 
K 3 C 60 . The evidence presented for superconductivity is that 
the instantaneous conductivity following photoexcitation in 
the normal state, (j(co’j), has features characteristic of the 
equilibrium (j(cd) measured at T C To obtain the in¬ 
stantaneous conductivity it is assumed that (7(cD;t) = (7(cd) + 
E((D;t), which is based on regarding E{(0\t) and 5a{(0\t) as 
equivalent. However, we have shown that these response func¬ 
tions are not equivalent, and differ most strongly when carrier 
momentum relaxation rates become long. Superconductivity, 
in which the condensate momentum relaxation time diverges, 
is the most extreme example of the inequivalence of E{co;t) 
and 5a{co;t). 

We considered two perspectives in attempting to model the 
response function E{co;t) appropriate to photoinduced tran¬ 
sient superconductivity. For a two-fiuid model with a tran¬ 
sient, fully coherent superfiuid component, we predict strong 
oscillations along the frequency axis of E{co;t), with a pe¬ 
riod inversely related to the sampling time. In the second per¬ 
spective, the pump induces a partially coherent state, with an 
enhanced, but still finite, momentum relaxation time. In this 
model, deviations from a Drude spectrum are again expected, 
though damped by the limited coherence time. In either case, 
the instantaneous conductivity G{(0\t) is unobservable or ill- 
defined in the most physically interesting regimes, while time- 
resolved THz-TDS measures E{co;t) directly. To advance re¬ 
search on photoinduced superconductivity and other collec¬ 
tive states, we believe it important to distinguish them. 
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Appendix A: Classical derivation of the nonequilibrium 
response 


To extend the classical Drude model to photoexcited mate¬ 
rials, we let both the carrier density n and the damping rate 
7 depend on time in the usual equation of motion. The cur¬ 
rent is then related to the field through the linear, first-order 
differential equation 


dJ / \ dn\ 
dt \ n dt) 


J = 



(Al) 


The extra damping term follows from the chain rule with 
J = nev, and causes the current to decay more rapidly when 
the carrier density decreases, as expected. Conversely, the 
damping term decreases when the carrier density increases, 
because our model incorrectly assumes that all carriers move 
with the same velocity.— Others have addressed this problem 
by expressing the current in terms of a distribution function, 
but their results reproduce Eq. (lAlb in the usual case of carrier 
decay.i2iii 

To solve Eq. dSB, we introduce the integrating factor 




njtj) 

n{t) 


exp[r(f ,-,0 {t-ti)], 


(A2) 


r(fi,?2) 


Ilfdt'rit') 

h-h 


(A3) 


the average damping rate over the interval (^ 1 ,^ 2 )- Assuming 
^ > 0 over we multiply Eq. (lAlb by (IA2I) and integrate 
to get 


J{t) —J{ti) = f dt'Qxp\—'Y{t',t){t — 

m Jti 

(A4) 

When n and 7 are constant, Eq. (IA4b gives the conventional 
Drude response. When n and 7 vary with time, the current 
J{t) includes contributions from impulses at earlier times t', 
exponentially weighted by the average damping rate experi¬ 
enced over its history. An impulse with carrier density n{t') 
will decay by a factor n{t)/n{t') before it contributes to the 
current/(r), so only the overall factor n{t) appears outside the 
integral. 

Referring now to Eq. •Ell, we let n = and 7 = 7 eq at 
equilibrium, and use Eq. (IA4I) to find the current change 5J{t) 
following photoexcitation aXt = 0: 


Z(T,r) = 0 (r) 0 (r - — {n{t)exp[-y{t - T,t)T\ 

m 

- neqexp(- 7 eqT)} . (A5) 

Eor the specific case of a photoinduced carrier density n{t) = 
^eq + with a constant Drude scattering rate 7 , 

Z(T,r) = 0 (r — T)S{t) exp(—rr) 0 (T)exp(— 7 t). 

m 

(A 6 ) 

The second factor of 0(r) in Eq. (IA 6 b is missing in Eq. (31) 
of Ref. [H but this appears to be a typographical error, since 
the factor is necessary to obtain their Eq. (32). Otherwise, the 
expressions are equivalent. 
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